1174 L.I. LAPIDUS 


cipally to the fact that Klein’s analysis (which _ OR. H. Dalitz, Proc. Phys. Soc. (London) A65, 175 

underestimates the role of states with L > 2) (1952). L. Wolfenstein and J. Ashkin, Phys. Rev. 85, 

agrees badly with the experimental data on n-p 947 (1952). 

scattering in the angular region close to 180°. 7], M. Shmushkevich, Thesis, Leningrad Physico- 
The author is thankful to V. P. Ddzelepov, Iu. M. ‘Technical Institute (1953). 

Kazarinov, and B. M. Golovin for profitable discus- 8 J. V. Lepore, Phys. Rev. 79, 137 (1950). 

sion and interest in the work. 9V. B. Berestetskii and I. Ia. Pomeranchuk, Dokl. 

Akad. Nauk SSSR 81, 1019 (1951). 

1}, Ia. Pomeranchuk, Dokl. Akad. Nauk SSSR 78, 249 10S. Courteney Wright, Phys. Rev. 99, 996 (1955). 
(1951). J. Exptl. Theoret. Phys. (U.S.S.R.) 21, 1119 11, Bakke and B. Stech, Z.Physik 144, 219 (1956). 
(1951). 121... Wolfenstein, Phys. Rev. 96, 1654 (1955). 

2G. Chew, Phys. Rev. 84, 710 (1951). 13 R. Oehme, Phys. Rev. 98, 147 216 (1955). 

3], M. Shmushkevich, Report to the Academy of Sci- 14, Klein, Nuovo cimento 2, 38 (1955). 


ences, U.S.S.R. (1953). 

4v. P. Dzhelepov, Thesis, Institute for Nuclear Prob- 
lems of the Academy of Sciences, U.S.S.R. (1954). 

5 Dzhelepov, Kazarinov and Fliagin, Dokl. Akad. Nauk =‘ Translated by E. S. Troubetskoy 
SSSR 100, 655 (1955). 284 


SOVIET PHYSICS JETP VOLUME 5, NUMBER 6 DECEMBER 15, 1957 


On the Magnetic Properties of Superconductors 
of the Second Group 


A. A. ABRIKOSOV 
Institute of Physical Problems, Academy of Sciences, U.S.S.R. 
(Submitted to JETP editor November 15, 1956) 
J. Exptl. Theoret. Phys. (U.S.S.R.) 32, 1442-1452 (June, 1957) 


A study is made of the magnetic properties of bulk superconductors for which the param- 


eter % of the Ginzburg-Landau theory is greater than 1/\/2 (superconductors of the second 
group). The results explain some of the experimental data on the behavior of superconduc- 
tive alloys in a magnetic field. 


Te AUTHOR’ has already noted that the quasi- Ginzburg and Landau considered only the case in 
microscopic Ginzburg-Landau theory” of super- which x «1/j2. 


conductivity leads to the conclusion that there ex- It has been shown by Zavaritskii,* however, that 
ist two groups of superconductors. For the first of _ the properties of pure metal thin films condensed at 
these groups, the parameter » entering into the liquid-helium temperatures are not described by 
Ginzburg-Landau theory is less than 1/V/2, and for such a theory. Zavaritskii and the present author 
the second group it is greater than 1//2. have therefore suggested that such films correspond 
This parameter x determines to a great extent to x > 1/2, and that superconductors can thus be 
the surface energy at the normal—superconducting ~— divided into two groups. The critical field for su- 
interface. It has already been mentioned,” in par- perconductors of the second group has already been 
ticular, that the calculated surface energy of a su- calculated’ as a function of the film thickness. The 
perconductor with x >1/\/2 is negative. Thus agreement obtained with Zavaritskii’s experimental 
superconductors of the second group should have data was not bad. 
properties very different from those of the first In the present work, a more detailed investiga- 
group. tion of the magnetic properties of bulk superconduc- 
For pure metal, x is found to be small. For in- tors of the second group (a cylinder in a longitudi- 


stance for mercury, x = 0.16. In view of this, nal field) is undertaken. The results obtained show 
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that x > 1/\/2 for a large number of superconduct- 
ing alloys whose magnetic properties had not pre- 
viously been well understood.* 


1, TRANSITION TO THE NORMAL STATE 


It was shown by Ginzburg and Landau? that if 
x > 1//2, the superconductivity is maintained at 
fields greater than H,,, at which equilibrium could 
exist between the normal and superconducting 
states. At fields higher than H,,,, a state with 
= 0 is unstable, and superconducting sections 
with & 40 may arise. It was shown’ that this in- 
stability continues to some value H,,, which for a 
bulk superconductor is x/2H,,,- At this value of 
the field the superconductor undergoes transition to 
the normal state by means of a second-order phase 
transition. In this section we shall investigate the 
properties of a superconductor in the neighborhood 
of the transition point, that is 

voy V 2A (1) 

The Ginzburg-Landau equations (in units previ- 
ously used by the author, and to be used hence- 
forth),t can be written : 


(<= st Ay = FPP, (2) 


—curl curlA =| PA + = (Y*VP—¥VY"). (3) 


We shall assume that the superconductor fills all 
space, and that the external field H, is directed 
along the z axis. Let the potential A be directed 
along the y axis. 

Close to the transition point | P|? <1, so that 
in the first approximation we may neglect the influ- 
ence of W on the field. We then obtain 

H—=const=H,, A= Hox. (4) 
We note that the point x = 0 may be located any- 
where in space. Let us now go on to Eq. (2). In- 
serting (4) into the equation for Y, neglecting the 
term | ¥ |? ¥ and considering ¥ to be a function of x 
only in the first approximation, we obtain the oscil- 
lator type of equation 


aE / dx? — x? (1 — Hox?) ¥ (x) = 0, (5) 


* The suggestion that x may be greater than W/2 for 
alloys was first made by L. D. Landau. 


{We note that in these units H,,, = 1/2. 
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* This equation has solutions 


when H, = x/(2n + 1). In particular, the largest 
value H, = x corresponds to 


as found previously. 


WY = exp {— 2x? / 2}, (6) 
In addition to such solutions, Eq. (2) is satisfied 
by the function 


exp fiky — = € — =) } (7) 


y2 


Since the point x = 0 has no special properties, the 
conditions will be exactly the same at all points of 
space, and it is natural to pick W in the form 


©O 
Y= Dy Creihy bn (x), 
n=-—-0O 


(8) 


with arbitrary coefficients / and C,. 

This form of V is a solution of the linear equa- 
tion, and refers in facts to 1, = x. In order to find 
a solution for H, < x, let us take into account the 
nonlinear terms in the equations. We find the ap- 
proximation of interest by inserting A from Kq. (4) 
with H, =, and W from Eq. (8) into Eq. (3). This 


leads to 


OPA 
Ox? 
= DEnCretir—mry [ice — CE") a (x) Ym (2), 
(9) 
= Beg = Be DEnCinei—O%9 (1 — m) Ya (2) Yo 2). 


It can be shown easily that 


H =0A/0x = H, —|¥ 2/2; (10) 


A = Hox —5.\| 8 Pdr, (11) 
is a solution of this equation, where V is given by 
Eq. (8), and H, is some still undetermined constant. 
It will be shown later that this constant is equal to 
the external field strength. 

Let us now consider Eq. (2) in the following ap- 
proximation. We insert into it Eq. (11), and take 
account of the term | ¥|*. We shall assume in this 
approximation that the function is still of the form 


of Eq. (8), but we shall add small terms Wx) to 
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all the C,y,,(x). The equations for these small terms are of the form 


PARGG) 
kn Bg 4 aryl) 4) nk 
(F— x] Yn — = a ata = 2 (x — Hy) »(x— 2) Cath (x) + 





(12) 


x 


+ 3) CpCnCnp4m{[ 2 — 5 (n— 25) | apm (4) | a2") Yn (2?) dx? — Bo (2) Yon 2) Yap tm (Oh 


p,m 





Thus yo? is a solution of an inhomogeneous equation. But in order for such a solution to exist, it is 
necessary that the right side of the equation be orthogonal to the solution of the corresponding homogene- 
ous equation. The solution of the homogeneous equation is just the function W,(x). Performing the neces- 
sary integration, we obtain the following condition: 


C,= 0. (13) 





(aT 1) > Co Caw Exp \— sa [(p —n)? + (m— ny} r : ‘a " 
p,m 


Let us multiply this equation by C* and sum over n. It is easily vernfied that we then obtain 





x— Hy 


% 





, 1 aes 
Pe + = a 1)(¥ z— 0, (14) 
where ¥ is the function of Eq. (8). Thus H, enters into W in the form of a coefficient. 


The quantity H, can be expressed in terms of the magnetic induction. According to (10) and (14), we 
have 


B=H=H,— | YP /2x = Hy—(x— Hy) / (22 — 1)8, (15) 
where 
B= Pw (FE ve 


is independent of H,. 
Let us now calculate the free energy. According to Ginzburg and Landau,? the free energy is 





Feopy—F 
SH so \ av \— | we + web 4 He 4 +(—! we —a¥) (Ye —AY*)}. (17) 
He (an “ * 
If we now make use of (2) and (3) and neglect the surface integral, which is of no importance to us, we 
obtain 
Fi.,—F 4 
= a3 \ (#2 — ar a —)av. (18) 
Hom | At , 
The free energy per unit volume is proportional to We note that (13) can be written in the form 
the average of the integrand. Denoting this aver- 
age by F, and writing it in terms of B, we obtain AB / OC%, 2G) (13') 
F 1 B (x — B)? (19) becoming equivalent to the condition that B is a 
13 ae 1 + (2%2—1)3° minimum. It is seen from (16) that 8 cannot be less 
: than unity. 
For a given & this quantity will be minimum when The condition that be periodic in x is the fol- 


6 has its lowest value. lowing simple recursion relation for the C,: 
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(20) 


Cine Cs 


The general analysis of the problem of choosing 
the coefficients is extremely complicated, but there 
is no reason to suppose that the actual value of N 
is a very large number. It would seem that the low- 
est energy corresponds to the simplest case* \ = 1. 
In this case B is given by 


Fk _ Ff 4\\ 
Pe (Ser | dye" }) 


The minimum value of B is obtained when k=xvV/2z, 


and is equal to 
y = > e—nn)? nee x (O, 1) = 1.18. (21) 


Then ¥ is 


yY =_ Ce-# 1294 [V 2a x1 (x -|- iy); 1]. (22) 


It can be shown from the properties of the function 
3, that rotation by 90° multiplies V by the phase 


factor e'~ ~*, but does not change it otherwise. 


*In addition to this case, which will be seen later to 
lead to a square lattice, we considered the example 
of a triangular lattice. The side of the triangle is 


(1/x) V 8x/3V 3. By numerical integration, however, 
we were able to show that this lattice corresponds to a 


larger value of 8, namely f = 1.32 (see also the next 
section). 


1177 


Thus ||? has the symmetry of a square lattice. 
When x = y = (1/x)/7/2, the function ¥ vanishes. 
According to (10), H attains its maximum value H, 
at these points. 

We shall now show that H, is equal to the field 
strength. As is well known, the field strength is 


42 OF /OB (in the usual units). This leads to a re- 
lation between H and B which is identical with Eq. 





(15). In order to perform the comparison with ex- 
periment (this will be done in Sec. 3) more conven- 
iently, the latter can be rewritten in terms of the 


function M(H,), 


— 42M | (Hep — Ho) = 1/1.18 (22 — 1), (23) 
where M and H can be taken in the usual units. We 
note that the magnetic moment is proportional to 
the first power of 1, — H,g. 


Equation (23) is the macroscopic characteriza- 
tion of the behavior of a bulk sample in an external 
field. As for the microscopic structure of the field 
penetration, collecting the results described above 
we obtain the following picture. In a plane perpen- 
dicular to the field direction, the structure is per- 
iodic with the symmetry of a square lattice. The 
field varies between its maximum value H, and its 
minimum value given by 


H min = Ho — (Hog — Ho) V2 [(22—1). (24) 
When the applied field strength //, is increased, the 
field between the maxima increases up to values 
equal to H =H, = %. 

The field and current distribution is shown in the 
figure. The lines correspond to | ¥ |? = const (nor- 
malized to | |4,,, = 1) or equivalently, to H 
= const. These are then the streamlines (their di- 
rections being given by the arrow). 


2, THE BEGINNING OF FIELD PENETRATION 
INTO A SUPERCONDUCTOR 


In the preceding section we concluded that in 
the neighborhood of the transition point to the nor- 
mal state, superconductors of the second group are 
in a special state which is neither normal nor super- 
conducting, and differs significantly from the inter- 
mediate state which may arise in noncylindrical 
bulk samples. Let us call this the mixed state. 
One naturally wonders how this state changes as 
the field is decreased, and where and how the tran- 
sition to superconductivity occurs. 

If we direct our attention to the experimental 
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B (H,) curves for alloys,* we see that this transi- 
tion can be both a first and a second order phase 
transition. It would seem (see Sec. 3) that the 
first order phase transition occurs when x is not 
very large, whereas the second order transition oc- 
curs when x >> 1. It is extremely difficult to ex- 
plain this on theoretical grounds, since this would 
require solving Eqs. (2) and (3) without any simpli- 
fying assumptions. 

Let us consider only the simplest case of a sec- 
ond order phase transition, which would seem, as 
has already been mentioned, to take place when x 
is not too small. 

Let us start with a determination of the transition 
point. It is natural to assume that when the field 
first begins to penetrate, the structure of this pene- 
tration is again of the lattice type, although with a 
very large period. At the lattice points the field 
will have its maximum value, and will practically 
vanish between them. Thus the field penetration 
pattern is that of separated filaments. This pat- 
tern reminds one of the distribution of superfluid- 
ity in helium II by the creation of vortex filaments, 
as proposed by Onsager® and Feynman.® The sim- 
ilarity is even greater when one analyzes the be- 
havior of the W function in the neighborhood of 
these filaments in greater detail. When one con- 
siders the situation at fields close to //,,, one 
sees easily that in the neighborhood of those points 
where W vanishes and the field is maximum, the 
phase of the ® function is* 6, where @ is the polar 
angle, and thus changes by 27 in going once around 
such a point. It is natural to suppose that the 
phase behaves in the same way in the neighborhood 
of separate filaments when the separation between 
them is large. As will be seen later, the W func- 


tion vanishes also in the center of a filament. This ~ 


behavior of the ¥ function, in particular of its 
phase, has been suggested by Feynman in applica- 
tion to vortices in helium II. 


When the filaments are sufficiently separated, 
they can be considered independent of each other, 
and each filament can be treated separately. In 
this case, the problem has cylindrical symmetry. 
Let us introduce the modulus and phase of the P 
function 


yr ae feis (25) 


* In the neighborhood of such a point with coordinates 
(x1, ¥1), ‘P is proportional to (x — 1) + ily — y,). 
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and consider the vector A directed perpendicular to 
the radius vector. We introduce the scalar quantity 
Q, which is equal to the absolute value of the vec- 


tor A—V/x. Equations (2) and (3) then become 


—a er f+ erat, 8 
alee (| = er. (27) 

The field H is given by 
H = ——-4 (rQ). (28) 


The statement of the problem requires that f> 1, 
“#20, G+O0asr-oo. Further, we have assumed 
that for a separate filament ¢ = 0. But this leads 
to Voy = dy/rd@ = 1/r. In other words, Q > 1/xr 
asr->0Q. According to the theory of Onsager and 
Feynman the superfluid velocity v, in helium II 
diminishes in the same way, with increasing dis- 
tance from the center of the filament. 

In addition, it is of course necessary to require 
that f be everywhere bounded. These conditions 
uniquely determine the solution of (26) and (27). In 
general, however, this solution can be obtained 
only by numerical integration, and this must be per- 
formed separately for each value of x. 

When x > 1, the solution is considerably simpli- 
fied. We note for this purpose that the distance 
over which Q should have a large variation is r ~ 1. 
As for f, it varies to a large extent at distances 
r~1/x. In view of this, one may assume that f 
has already attained the value f = 1 in Eq. (26), 
from which we immediately obtain 


Q = Ky (r)/x, (29) 
where K, is the Hankel function of imaginary argu- 
ment. Further, we note that Q = 1/xr whenr <1. 
Inserting this into the equation for f we obtain 


a [— & (r a = — | — fs—f, (30) 


At distances r > 1/x, the solution of this equation 
is f? = 1 — 1/x?r?, and when r <1/x, the solution 
is f = Cr, where C is a constant which can be ob- 
tained by numerical solution of (30). 

Let us now calculate e, the free energy per unit 
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length of a filament. To do this we must evaluate 


the integral 


c= 2s (H+! 


0 





<5) rar 


> | (31) 
= | (1 a f+) r— p2 <. (72Q2)| dr. 


If we note that the integrand decays exponentially 
when 7 > 1, and assume that In x > 1, it is easily 
seen that the important contribution to the integral 
will be from the interval 1/x «<r<1, and we then 
obtain 


c= (27/2) Inx. 


Numerical integration of (30) and (31) makes it pos- 
sible to obtain a correction to In x, so that we fi- 
nally obtain 
e = (27 / x?) (Inx + 0.081). (32) 

In order that the formation of filaments be ener- 
getically possible, the field strength //, must be 
such that F, — 2H,B be negative. If the filament 
density, i.e., the number of filaments crossing a 
unit area, be denoted by n, we have 
Flan, B=H=n | HdS =n 6 Adi. (33) 
In view of the fact that the contour of integration 
extends to r = oo, and that at these large distances 
Q =0, we conclude that A= Vy/x, which means 
Thus 


that its absolute value is 1/xr. 


B= 2nn/x. (34) 
It follows from this that the filaments start being 
formed when the field strength //, attains the value 
Hey = ex / An. (35) 

This formula is valid for arbitrary x on the assump- 
tion that the phase transition is of second order. 
When x > 1, we obtain 

Ag, = (1 / 2x) (Inx + 0.08). (36) 

It is interesting to note that at the center of the 

filament the field is about twice H,, in the given 


Indeed, according to Eq. (27) we have 


case. 
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\ Qfedr = (Inx — 0,18) /x. 


0 


H (0) = (37) 


From (35) it is seen that ,, is independent of n, 
and that an arbitrary number of filaments may be 
formed in this approximation for H, > H,,- In order 
to find the function n(H,), and therefore B(H,), we 
must take account of the interaction between the 
filaments. 

This is easily done if x > 1. Indeed, in this 
case the region in which f differs from unity (r 
~ 1/x) is small compared with the distances at 
which the interaction between the filaments takes 
place (r ~ 1). In view of this, the filament cen- 
ters will play the role of singularities in the solu- 
tion, and their detailed structure will not be of im- 
It is not difficult to see that the equa- 
tion for the field together with the conditions at 
the centers of the filaments can be written in the 
form 


portance. 


AH — H = — (2n/%) 8 (r —tm), (38) 


where the r,, are the coordinates of the centers of 


the filaments. The solution of this equation is 


1 
H =—- DiKo(\t —tm). (39) 
The energy can be written in the form 
F, = A? + (VH)? = H (H — ABH) 
(40) 


= 2S Ko(\tml); 


where the r,, are calculated from one of the filament 
centers. The term corresponding to r,, = 0 should 
be replaced by Eq. (32). The quantity n and the 
periodicity of the lattice, which enter into (40), can 
be expressed in terms of the magnetic induction 5 


by Eq. (34). For a square lattice we obtain 


Bet bY KY 2EE™) (an 


% 
Fy — 9 
[4-m2>1 


It is, however, possible that the lattice is actu- 
ally not a square one or that its symmetry changes 
as the field is varied. Therefore the problem 
should be dealt with more carefully. The most 
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probable structure which one can assume instead 
of the square one is the triangular one. For such a 


lattice one — obtains 


¥B An (2 + m? + Im)’ 
f,= a On € + > Kt OV 8yR xB ast 
24+ m?+Ilm>1 (42) 


By comparing the asymptotic form of this expression 
for small B with asymptotic form of (41), 


xB. 
20 


eV eae) ex (-V 5): 


pe 
B 4B 
From pet 42V F( 


one sees clearly that for sufficiently small B there 


Eines 





(43) 


should be a transition to the triangular modification. 
The value of H, at the transition point is deter- 
mined by setting the free energy for a given field 
strength equal to fF, — 2H,B, and is 
Hy = Hey + 0.0394 /x. (44) 
This transition is clearly a first order phase tran- 
sition. The induction undergoes a jump from B, 
= 0.286/% to B, = 0.294/%. At fields greater than 
H;, the square lattice is more favorable, which jus- 
tifies to some degree the assumption made in the 
previous section. 

We note that the transition takes place in the im- 
mediate neighborhood of //,,, and that the discon- 
tinuity in the induction is so small (about 3%) that 
it would be extremely difficult to observe it on a 
B(H,) curve. 

It is not difficult to see that when // = H,,, the 
B(H,) curve has a vertical tangent. Indeed, when 
H = H,,, the 1,-dependence of B is determined pri- 
marily by the term H, ~ expl| (42 /\/3xB)? |, and 
dH,/dB > 0 as B > 0. In the other limiting case, 
when x5 > 1, the various terms of the sum of Eq. 
(41) have arguments which lie close to each other, 
so that the sum can be replaced by an integral. 
This gives F, ~ B’, andH, =~ B. Thus as UH, is in- 
creased, the curve asymptotically approaches the 
straight line 5 =H,. This goes on until the dis- 
tance between the filaments becomes of the order of 
1/x, i.e., until H, ~ x. 

The dependence of xB on (H, — H,,) is given by 


the universal function 


se) exp(—V 3), 
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% (L719 — Hex) = —y oo (X1,m) ONE, mi (X1,m)]; 
X1,m 


+ m? + lm) /V3xB 


x (Hp — Ha) < 0.0394, 
Xim = V 2a (22 + m2) /«B 


x (Hy — He) > 0.0394. 


where X1,m = V 4x (1? 


(45)’ 


We note that (0B/OH,)7 is always positive, which 
indicates that the situation discussed is stable. ’ 

For very large values of x Eq. (35) can be used 
to find H,,. To do this one must solve (26) and 
(27) numerically, after which © is calculated from 
(31). On the other hand, finding B(H,) in this case 
requires solving a separate and very complex prob- 
lem. 

To conclude this section, let us turn our atten- 
tion to the intermediate state in noncylindrical 
speciments (or cylinders in nonparallel fields). It 
is most interesting that no intermediate state arises 
in specimens for which x > 1. One obtains instead 
a uniform pattern of a mixed state (in ellipsoidal 
samples). When the transition in the field H,, is of : 
the first order, however, an intermediate state does 
arise. 


3. COMPARISON OF THEORY AND 
EXPERIMENT 


No detailed data exist currently on the magnetic 
properties of superconducting alloys. The only 
work which is more or less suitable for comparison 
with the theory is that of Shubnikov, Khotkevich, 
performed in 1937. Let 
us compare the theoretical results with the data of 
this work. 


Shepelev, and Riabinin, * 


If we are to speak of general qualitative charac- 
teristics, we must first note that starting at some 
alloy concentration, the B(//,) curve indicates a re- 
gion of partial field penetration; this corresponds 
to the conclusion reached above on the existence 
of a mixed state. This region is bounded by the 
two field strengths H,, and H,,, and although the 
transition from the superconducting state to the 
mixed one can be either of first or of second order, 
the transition to the normal state must definitely be 
of first order. In the neighborhood of H,,, the 
B(H,) curve is linear. All this is in agreement with 
the results of Sec. 1. Further, if the alloy concen- 
tration is sufficiently high, the transition at H,, 


seems also to be of second order. According to the 
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Alloy ; H H tan O tan O Hey Hes 
composition rons - - (theory) | (exp.) Gis) (exp.) 
Pure Pb 4h | 540 | 0.22-+0. on puers — [a| | swoforroxl = | ff | 
SP aes ae Na ae ie tha Ss aa em ea ee oe as ee Se ee 
4.22} 704 | 596 1.18 2.12 2.2 ss 540 
Pb+2.5% Tl | 4/99 | 14100 | 846 1.3 1.24 1.2 i 740 
4.22} 1010 | 564 1.8 0.38 0.38 - 430 
Poh +5% Tl | 4°89 | 4730 | 354 9 0.27 0.27 600 
4.22] 2470 | 534 4.07 0.054 | 0.052 150 230 
Pb + 15% Tl | 4°76 | 3760 | 762 4.9 0036 | 0.036 210 310 
4.22, | 2840 | 575 4.94 0.036 | 0.033 150 170 
Pb + 30% Tl | 4°79 | 4460 | 756 5c0 0.025 | 0.023 190 240 
4.22. | 2270 | 370 6.14 0.023 | 0.02 90 140 
Pb + 50% Tl | 93° | 470 | 520 8.2 0.013 | 0.044 120 150 
4.22,| 983 | 604 1.6 0.54 0.45 ae 530 
Pbh-+2% In | 4'95 | 4520 | 342 1.8 0.37 0.3 a 730 
4.22 | 2400 | 563 4. 26 0.049 | 0.045 170 220 
Pb + 8% In 1.75 | 3780 | 842 4.7 0.044 | 0.043 240 310 


theory, at this point the B(H,) curve has a vertical 
tangent, and then begins gradually to approach the 
straight line B = //,. This conclusion also more or 
less agrees with the experimental curves. 

Going on to a quantitative comparison, we note 
that the two parameters H,,, and x which enter into 

the theory must be determined from the experimental 
curves. This can be done by finding the upper crit- 
ical field H,, and the area under the curve of the 
magnetic moment B(H,) — H,, where H, is the field 
strength. In fact from the definition of H,,,, we 
have 


4H, = (Hy —B (A) dite 


Ce 9 Q 


HT ¢9 
= \ [H.— B(H)] dH. 


0 


(46) 


Of course in order that /],,, defined in this way 
have meaning, it is necessary that the experimental 
B(H,) curve be that for equilibrium. When we know 


Hem» We can find x from the equation H,, = xV 2H, m- 


The values of H.,, 
can be used to compare the magnetic moment in the 
neighborhood of the upper critical field H,, with ex- 
periment. According to Sec. 1, the dependence of 
M on H, is given by Eq. (23). 

We can compare with experiment, in addition to 
M, Eq. (36) for the first critical field. 


The results of such comparison are given in the 


and » obtained in this way 


Table.* One should bear in mind that the accuracy 
with which the slope of the magnetic moment curve 
is determined in the neighborhood of H,, (let us de- 
note the ratio in Eq.(23) by tana) and H,, 


curves of Shubnikov, et al.4 


from the 
is no greater than 15%. 
It is seen from the table that the data for tana are 
in very good agreement with the theory. The values 
of tana vary from specimen to specimen by a fac- 
tor as high as 150, and the theory always gives the 
correct value within the above-mentioned limits. 
The situation with respect to H,, is somewhat 
worse. Although the theoretical values are always 
of the same order of magnitude as the experimental 
ones, they lie too low in all cases. It should be 
borne in mind, however, that the values of x dealt 
with here are not very large compared to unity, and 
this can lead to a systematic error in determining 
H., from (35), an equation which is valid only when 
x >> 1. It is not difficult to see in which direction 
the theoretical value must be corrected in order to 
account for the fact that x is not very large. This 
can be done by making the natural assumption that 
Ho: =Hem when x = 1//2. Instead of this, inser- 
tion of x = 1//2 into Eq. (36) leads to a negative 


value of H,,. Thus for values of x which are not 


* It is significant here that the theory remains appli- 
cable for temperatures far from T,. As fields close to 
Hea» this conclusion is strictly correct, since P is small 
in this region. Recently, however, Ginzburg® has pre- 
sented arguments which tend to show that the theory is 
applicable also for arbitrary fields and temperatures. 


1182 


very large, the actual theoretical value of H,, 
should be greater than that given by (36). In prin- 
ciple, it is also possible to perform a comparison 
between experiment and Eq. (45) for B(H,) in the 
neighborhood of H,,. This equation, however, is 
in principle valid only for x > 1, so that only a 
qualitative comparison is possible; this has already 
been done at the beginning of this section. Further- 
more, the experimental curves in this region are 
measured with very low accuracy, so that a quanti- 
tative comparison is not really very meaningful. We 
note here that according to the experimental data, 
when x is close to 1/2 the transition to the mixed 
state would seem to be a first order phase transi- 
tion. In this case the concept of isolated filaments 
does not correspond at all to the mixed state at the 
transition point. | 

In order to investigate the dependence of #,,, 
and » on the concentration, let us compare their 
values at a fixed temperature of 4.2°. For compari- 
son, the table gives the data also for pure lead at 
the same temperature, taken from Lock.” It is seen 
from the table that H,,, hardly depends on the con- 
centration (except for the 50% Pb-Tl alloy). The 
variation is no greater than 11% and has no appar- 
ent regularity. However, x increases regularly with 
‘increasing concentration, and varies for Pb-T] al- 


loys from 1.18 at 2.5% TI to 6.14 at 50% TI. 


Some years ago Pippard?® made the interesting 
suggestion that there exists in superconductors a 
certain characteristic length (the coherence length) 
which determines the “range of order” of the super- 
conducting electrons. This length characterizes 
the minimum distance at which the number of su- 
perconducting electrons may change significantly. 
It may, for instance, refer to the thickness of the 
transition layer between the normal and supercon- 
ducting phases. 

According to the theory of the previous sections, 
the quantity 6/x should play the role of such a 
length. This length determines the thickness of 
the transition layer in the intermediate state.” If 
we bear in mind that x = (\/2e/hc)H, 5, then it is 
clear that the characteristic length is proportional 
to 1//x. According to Pippard, this length should 
depend strongly on the impurity concentration. At 
high concentrations it should be of the order of the 
mean free path of an electron in the metal. It is 
seen from the table that this concept is not in dis- 
agreement with the experimental data. 
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The results described show that the quasi-micro- 
scopic theory is able to describe many of the regu- 
lar properties of superconducting alloys. Thus 
these regular properties arise primarily not because 
macroscopic regions with relatively large values of 
[, and H.m arise in the alloy, as has been thought 
previously, but from the destruction of electron 
order at distances small compared with the pene- 
tration thickness 6. | 

It is true, of course, that macroscopic inhomo- 
geneities may be of importance. Without even go- 
ing into the fact that a strong inhomogeneity can 
change the overall pattern, we note that such a phe- 
nomenon as the remanent magnetic moment cannot 
be described by the present theory and is more 
probably a result of such inhomogeneities. The 
structure of the field penetration itself should be 
extremely sensitive to inhomogeneities, so that in 
real substances it would hardly be in the form of a 
regular lattice. 

In conclusion [| take this opportunity to express 
my gratitude to Academician L. D. Landau and N. 
V. Zavaritskii for valuable comments and discus- 
sion of the work. 
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